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ABSTRACT 

The a Cen AB system is an attractive one for radial velocity observations to detect poten- 
tial exoplanets. The high metallicity of both a Cen A and B suggest that they could have 
possessed circumstellar discs capable of forming planets. As the closest star system to the 
Sun, with well over a century of accurate astrometric measurements (and a Cen B exhibiting 
low chromospheric activity) high precision surveys of a Cen B's potential exoplanetary sys- 
tem are possible with relatively cheap instrumentation. Authors studying habitability in this 
system typically adopt habitable zones (HZs) based on global radiative balance models that 
neglect the radiative perturbations of a Cen A. 

We investigate the habitability of planets around a Cen B using ID latitudinal energy 
balance models (LEBMs), which fully incorporate the presence of a Cen A as a means of 
astronomically forcing terrestrial planet climates. We find that the extent of the HZ is relatively 
unchanged by the presence of a Cen A, but there are variations in fractional habitability 
for planets orbiting at the boundaries of the zone due to a Cen A, even in the case of zero 
eccentricity. Temperature oscillations of a few K can be observed at all planetary orbits, the 
strength of which varies with the planet's ocean fraction and obliquity. 
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1 INTRODUCTION 

There are currently 755 known extrasolar planets, spanning a large 
parameter space in orbital elements and in planetary propertieo 
The question of planet formation in binary systems is an impor- 
tant one, as a large fract ion of solar type stars are b om in binary 
or multiple star systems touguenn ov & Mavoi)ll99ll) . The major- 
ity of these systems are wide S-type binaries, where the secondary 
orbits at distances of over 100 au from the planetary system. At the 
time of writi ng, around 15% of planets h ave been detected in binary 
star systems besidera & Barbieril ( l2007h . In S-type systems, planet 
formation modes should be similar to those in single star systems, 
with the secondary's influence being roughly negligible (unless its 
orbit is sufficiently eccentric). 

However, there are a handful of systems with detected planets 
where the se condary orbits at m uch closer se parations (^ 20 au ) 
- 7 Cephei b ('Hatzes et al."2003), HD41004b I zucker et al.ll2004h . 
GJ86b IqucIoz et al. 2000), etc. The a Centauri system has similar 
orbital architecture, although as yet it does not host detections of 
exoplanets. 

Despite this, the a Cen system is an attractive one for studying 
possible planet formation in binary systems. At a distance of 1.33 
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pc, it is the nearest star system to our Sun. It is composed of a hier- 
archical triple system, with the central binary system (a Cen AB) 
orbited by the M dwarf Proxima Cen at suffic iently large distance to 
be negligible jWertheimer & LaughUij2006h . a Cen A is type G2 V 
with mass Ma = 1.105 ± O.OOTM©, and a Cen B is type KIV 
with mass Mb ~ 0.934 ± O.OOTAf©. It is particularly amenable 
to rel atively cheap radial velocity (RV) campaigns ( Gued es et al.l 
|2008|) . especially if the goal is to detect Earth-mass planets within 
the habitable zone of a solar-type star. 



As iGuedes et alj J^OOt) describe, there are several reasons 
why this is the case: firstly, both a Cen A and B are high metallic- 
ity stars, which would promote the existence of circumstell ar discs 
with a high fraction of solid materials at early times (Wv att et"al] 
l2007h . as well as deeper spectral lines for improved RV precision. 
Q Cen B is particularly quiet in terms of chromospheric activity 
and acoustic p-wave oscillations, with a relatively strong potential 
RV signal due to its lower mass. The binary is inclined by only 1 1 
degrees with respect to the line of sight, i.e. its inclination angle 
from face-on is i = 79°. This ensures that sini ~ 1, and that the 
recovered planet mass from RV observations will be close to the 
true mass, provided that the planets form and remain in the same 
plane. Finally, its position in the sky (—60° declination) affords as- 
tronomers in the southern hemisphere the opportunity to observe 
the system for most of the year. If a Cen B does host terrestrial 
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planets, they should be readily detectable even with a 1 metre tele- 
scope with high-resolution spectrograph instrumentation. 

Several numerical studies have shown that both stars are ca- 
pable of forming terrestrial planets despite the perturbing influence 
of the binary companion, which appears to play a role analogous 
to the gas giants within our solar system. The planetesimal discs 
appear to be stable within 3 AU of their parent stars, provided the 
inclin ation of the disc relative to the binary plane is less than 60 de- 
grees dWiegert & Holmanlll99llOuintana et alj2002l l2007'). How- 
ever, other studies have shown that the later stages of accretion to 
produce lunar mass objects is reduced in efficiency due to orbital 
rephasing by the binary companion. This inhi bits collisional growt h 
around a Cen A to regions w ithin 0.75 au {Thebault et al.ll2008h . 
and within 0.5 au of a Cen B jThebault et al Il2009h . 

It therefore appears to be the case t hat gas giant fo rmation 
is suppressed relative to our solar system jXie et allllOlOl) . a pre- 
diction consistent with the absence of detections from previous ra- 
dial velocity surveys, which suggest upper limits between 0.5 and 
3 Mju p. If the disc can produce a sufficient quota of lunar mass 
bodies. lOuedes et aT, 12008) show that Earth mass planets can form 
in the habitable zone of a Cen B, with maximum eccentricities of 
around ep = 0.3 (this result is also corroborated by Xie et al 2010). 

Once formed, Hamiltonian analysis indicates terrestrial plan- 
ets orbiting in a Cen B's habitable zone appear to be dynamically 
stable under perturbations from a Cen A, provided that < 0.3 
and the inclin ation of the planet's orbit relative to the binary plane, 
ip < 35° jMichtchenko & Porto de Melldl2009l) . Equally, these 
authors also show that planets with inclinations larger than this 
value are expected to experience strong instability due to the Lidov- 
Kozai resonance resulting in eccentricity-inclination coupling. 

These studies have generally assumed a habitable zone (HZ) 
around a Cen B in the semi major axis range 0.5 < a < .9 au. 
This range is based on calculations by Kastin g^t alj l ll993h . who 
calculate the HZ using a global radiative balance (GRBM) model, 
assuming Earth mass planets with similar atmospheric composi- 
tion (N2/H2O/CO2). The inner edge of the HZ is governed by loss 
of water via photolysis and subsequent hydrogen escape, and the 
outer edge is determined by formation of CO2 clouds, which in- 
crease the planet's albedo. In defining the habitable zone in this 
fashion, the perturbing influence of a Cen A is neglected. Given 
the semimajor axis of the orbit, this would appear to be an appro- 
priate approximation. If main sequence relations for the luminosity 
of each object are assumed, the insolation experienced by planets 
in the habitable zone of a Cen B due to a Cen A would be no more 
than a few percent of the total insolation of the a Cen AB system at 
the binary's periastron, and around one tenth of a percent at apas- 
tron. This insolation can be diminished further by eclipses of a Cen 
A by Q Cen B, the duration of which is estimated to be of order a 
few Earth days. 

However, we should also consider the results of more com- 
plex ID latitudinal energy ba lance models (LEBMs) such as 
those describ e d by [ Williams & Kastind 1 1997 ) an d subsequently 
ISpiegel et al.l j2008l ,2009) ; Dressing et alj ilOld) : ISpiegel et al.l 
( I2OIOI) . Rather than pursuing a simple global balance, the planet's 
temperature is allowed to vary as a function of latitude, A. The in- 
solation of the planet will also be a function of latitude and season, 
and the other key properties (infrared cooling rate and albedo) are 
also temperature dependent. 

A planet in global radiative balance is not in general in lo- 
cal radiative balance, and by extension habitability is not a discrete 
concept (i.e. either habitable or uninhabitable), but a continuous 
one, where a certain fraction of the planet's surface will be habit- 



able at any given time. In the LEBM, the evolution of the planet's 
temperature r(A) is described by a diffusion equation made non- 
linear by the addition of the heating and cooling terms, as well as an 
albedo which makes a rapid transition from low to high as temper- 
ature decreases past the freezing point of water. As a result, small 
changes in the properties of a planet can strongly affect the resultant 
climate. For example, changing the length of day in an Earth-type 
model can determine whether the planet can retain liquid water on 
the surface, or undergo an albedo feedback reaction which results 
in the "Snowball Earth" scenario dSpiegel et al1l2008l) . 

These models have been successful in establishing important 
aspe cts of astronomical forcings on climate, e.g. Milankovitch cy- 
cles ( ISpiegeIetall2010l) . variations due to orbital eccentricity, and 
potentially the effects of Kozai resonances or other orbital instabil- 
ities. It is this propensity for incorporating external forcing (such as 
the presence of a binary companion) which lends it towards study- 
ing habitability in systems such as a Cen AB. The sensitivity of 
the climate to such forcings suggest that the relatively small pertur- 
bation produced by the presence of a Cen A may have important 
consequences for the location of the habitable zone around a Cen 
B. 

In this work, we augment the ID LEBM of ISpiegel et al.l 

( |2008[) to include the presence of a binary companion, and perform 
a parameter study for planets orbiting a Cen B. In particular, we in- 
vestigate the currently defined habitab l e zone , to compare with the 
GRBM calculations of iKasting et al.l ( Il993l) . Section |2] describes 
our modified LEBM, and the initial conditions studied in the pa- 
rameter space survey. Section [3] displays the results of this study. 
In section |4] we discuss the implications of these results, and sum- 
marise the work in section[5] 



2 METHOD 

2.1 Latitudinal Equilibrium Balance Models 

At their core, LEBMs solve the following diffusion equation: 

where T — T{x,t) is the temperature at x = sin A, and A is the 
latitude (between —90° and 90°). This equation is evolved with 
the boundary condition ^ = at the poles. The (1 — x^) term 
is a geometric factor, arising from solving the diffusion equation in 
spherical geometry. 

C is the effective heat capacity of the atmosphere, D is a diffu- 
sion coefficient that determines the efficiency of heat redistribution 
across latitudes, S is the insolation flux, / is the IR cooling and A 
is the albedo. In the above equation, C,S,I and A are functions of 
X (either explicitly, as S is, or implicitly through T{x)). 

D is a constant, defined such that a planet at 1 au around a 
star of IMq, with diurnal period of 1 day will reproduce the av- 
erage temperature profile measured on Earth. Planets that rapidly 
rotate experie nce inhibited latitudinal heat transport, due to Corio- 
lis effects (see lFarrelllll990l) . In the model, we follow [Spiegel et al] 
( 2008) by scaling D according to: 



D = 5.394 X 10^ 



(2) 



where ujd is the rotational angular velocity of the planet, and LJd,(B 
is the rotational angular velocity of the Earth. This expression is 
probably too simple to describe the full effects of rotation, but in the 
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absence of a first-principle theory to describe that describes reduced 
transport we must make do. 

In this work, we solve the diffusion equation using a simple explicit 
forward time, centre space finite difference algorithm. A global 
timestep was adopted, with constraint 



St < 



2D{1 



(3) 



The parameters are diurnally averaged, i.e. a key assumption 
of the model is that the planets rotate sufficiently quickly relative 
to their orbita l period. We adopt the same expressions for these 
parameters as ISpieg el et al. I ( l2008h . who in turn used the work of 
IWilliams & Kastin g|( ll997h . 

The atmospheric heat capacity depends on what fraction of 
the planet's surface is ocean, f ocean, what fraction is land fiand ~ 
1.0 — f ocean, and what fraction of the ocean is frozen fice'. 



C — flandCland 4" f ocean ((1 fice^Cocean ~t~ ficeCice) • 

The heat capacities of land, ocean and ice covered areas are 
Ciand = 5.25 X 10" erg cm-^ 

Co cean — and 



Cice 



9.2Cland 
2Cland 



263 K < T < 273 K 
T < 263 K. 



The infrared cooling function is 

UT) = '^^^ 

^ ' l + 0.75rjfl(r)' 

where the optical depth of the atmosphere 

riR{T) ^ ^ ^ 



0.79 



^273K 
The albedo function is 

A{T) = 0.525 - 0.245 tanh 



r - 268 K 



5K 



(4) 

(5) 
(6) 

(7) 

(8) 

(9) 
(10) 



This produces a rapid shift from low albedo to high albedo as the 
temperature drops below the freezing point of water. It is this tran- 
sition that makes the outer habitable zone extremely sensitive to 
changes to various orbital and planetary parameters, and makes 
LEBMs an important tool in studying short-term temporal evolu- 
tion of planetary climates. The insolation flux 5* is a function of 
both season and latitude. At any instant, the bolometric flux re- 
ceived at a given latitude at an orbital distance r is 



S = qo cos Z 



(11) 



where go is the bolometric flux received from the star at a distance 
of 1 AU, and Z is the zenith angle: 



go = 1.36 X 10" 



M 



-1 -2 

ergs cm 



cos Z — fj. — sin A sin S + cos A cos S cos h. 



(12) 
(13) 



We have assumed main sequence scaling for the luminosity (AIq 
represents one solar mass). Given that both a Cen A and B are 
similar in mass (and spectral type) to the Sun, this is a reasonable 
first approximation, however, we should note the observational con- 
straints placed by Theve nin et al. U2003) . as we will in the Discus- 
sion. S is the solar declination, and h is the solar hour angle. The 
solar declination is calculated from the obliquity So as: 



sin 5 = — sin So cos{(j>p — (f)peri — (f'a), 



(14) 



where (/)p is the current orbital longitude of the planet, (jjperi is the 
longitude of periastron, and is the longitude of winter solstice, 
relative to the longitude of periastron. 

We must diurnally average the solar flux: 



(15) 



This means we must first integrate /i over the sunlit part of the 
day, i.e. h = [~H,+H], where H is the radian half-day length 
at a given latitude. Multiplying by the factor H /n (as H — vr if 
a latitude is illuminated for a full rotation) gives the total diurnal 
insolation as 

S = qo ( —] p. = — {H sin A sin S + cos A cos S sin H) . (16) 

\ TT J TT 



The radian half day length is calculated as 
cos H — — tan A tan S. 



(17) 



We calculat e hab itability indices in the same manner as 
ISpiegel etal] ( l2008l) . The habitability function 77 is: 



7?(A,t) = 







273 K < T{X,t) < 373 K 
otherwise. 



(18) 



We then average this over latitude to calculate the fraction of hab- 
itable surface at any timestep: 



(19) 



We will use this function to classify the planets we simulate in the 
following sections. 



2.2 Modifications to include tiie binary 

The addition of the second star requires us to repeat the insolation 
flux calculation, where we must now re-calculate the orbital lon- 
gitude, solar declination and radian half-day length of the planet 
relative to the secondary, and use the current the orbital distance 
from the planet to the secondary. 

We must also account for the possibility of transits. Given the 
orbital configuration, transits of the secondary (a Cen A) by the 
primary (a Cen B) should occur frequently in the simulation. By 
calculating the angle between the vector r2i between the two stars, 
and the vector rpi between the planet and the primary, it can be 
determined whether a transit is occurring at any given timestep 0. 
If the primary transits the secondary, the secondary flux is set to 
zero. 

We ensure that transits are resolved in time by adding a second 
timestep criterion, which ensures that at the planet's current orbital 
velocity, the duration of a transit will not be less than 10 timesteps. 
We neglect the distant companion Proxima Centauri. 



2.3 Initial Conditions 

Unless otherwise stated. Earthlike conditions are assumed. The di- 
urnal period is equal to the Earth's, and the obliquity is set to 23.5° . 

f ocean IS Set tO 0.7. 

The a Cen system is set up to have Mi = Mb ~ 0.934Mq, 



^ We also assume the stellar radii are governed by main sequence relations, 
see lPrialnilj i2OO0l) 
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M2 = Ma = I.IMq. The semimajor axis of the orbit is a = 23.4 
au, and the eccentricity e = 0.5179. a Cen A is placed at apastron 
at the beginning of all simulati ons. 

In line with the results o f iGuedes et alj ( I2OO8I) , we investigate 
planet eccentricities between = to Cp = 0.3, and semi major 
axes relative to ct Cen B in the conventionally established habitable 
zone 0.5 < a < 0.9. 

The simulations begin at the northern winter solstice, which 
is assumed to occur at an orbital longitude of 0°. In the case of 
eccentric orbits, this is also the longitude of periastror[f[ The plan- 
ets' initial temperature was 288 K globally (starting the simulation 
at higher temperatures did not significantly affect the result). The 
simulations were carried out for 1000 years (approximately 15 or- 
bits of Q Cen A around a Cen B). The first 200 years of simulation 
are ignored, allowing the system to settle down to a quasi-steady 
state in all cases. 
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Figure 1. The habitable zone around a Cen B. The x-axis shows the semi- 
major axis of the planet orbiting a Cen B, the y-axis shows the eccentricity 
of the planet's orbit. The resulting planets are classified into five categories, 
which are described in the text. 



3 RESULTS 

3.0. 1 The Habitable Zone 

Figure [T] shows the end result for all simulations carried out in this 
work. The various end states can be classified thus: 

(i) Habitable Planets - these planets are 100% habitable across 
the entire surface, i.e. f^{t) — 1 for all t. 

(ii) Hot Planets - these planets have temperatures above 373 K 
across all seasons, and are therefore uninhabitable {fj{t) = for all 
t). 

(iii) Snowball Planets - these planets have undergone a snow- 
ball transition to a state where the entire planet is frozen, and are 
therefore uninhabitable (f]{t) — for all t). 

(iv) Eccentric Transient Planets - these planets are partially hab- 
itable (0 < fj{t) < 1), but the habitability fraction oscillates ac- 
cording to the period of the planet's orbit around a Cen B. This 
oscillation is not present in circular orbits. 

(v) Binary Transient Planets - these planets are partially habit- 
able (0 < f]{t) < 1), but the habitability oscillates with the period 
of a Cen A. This oscillation is present in both circular and eccentric 
orbits. 

Some systems exhibit both types of transient behaviour Where this 
is the case, we classify systems according to which oscillation has 
the largest amplitude. 

If we instead chose to classify by maximum temperature, then 
most systems would be regarded as binary transients. Even plan- 
ets in the middle of the habitable zone will experience temperature 
fluctuations. The resulting temperature oscillations due to a Cen A 
are quite small - the maximum temperature generally fluctuates by 
at most 2K. 

Figure [2] shows the total fraction of the planet which is habitable, 
f]{t), time-averaged over the range [300, 400] years. The steep gra- 
dient at the outer edge of the HZ indicates the non-linear, highly 
sensitive nature of the snowball transition. For single-star models 
using parameters corresponding to Earth and the Sun, ~ 0.85. 
Around a Cen B, planets on circular orbits will attain similar val- 
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Figure 2. Habitability Fraction as a function of planet semi-major axis, time 
averaged in the date range [300, 400] years. Curves are plotted for four 
values of the planet's orbital eccentricity. 



ues at Bp = 0.86 au0. This plot would indicate that the inner HZ 
would begin nearer to 0.6 than 0.5 au, with the edge of the outer 
HZ agreeing with previous results. 



3.1 Habitable, Hot and Snowball Planets 

Planets in these three categories have stable values of rj throughout 
the season. The influence of a Cen A tends to raise the temperature 
by a few K at all eccentricities, with this phenomenon insensitive to 
any resonances between the longitude of periastron of the star and 
planet. 

Temperature profiles (time averaged over one orbit of the 
planet) take essentially the same shape, of the form 



^ Simulations were carried out where the longitude of periastron was var- 
ied. The results were not significantly affected, as the eccentricities studied 



are relatively low, see [Dressing 



cantly attecte' 



T{\) ^ A~ Bsin X 



(20) 



With the snowball planets taking the lowest value of B, and the hot 
planets taking the largest values of B. 
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Figure 3. Habitability fraction as a function of time for tlie case where 
ap = 0.675, Cp = 0.3. 

3.2 Eccentric Transience 

As was found bv lPressing et aljfeoioh . eccentric planets are gener- 
ally hotter compared to planets in circular orbits at the same semi- 
major axis (as received flux scales as (1 — e^)"^''^). As the eccen- 
tricity regime we explore is fairly modest, we do not see planets 
entering and leaving the snowball state purely because of high Cp. 
The planets do not spend long time intervals at apastron (relatively 
speaking), and therefore would not be able to freeze except if the 
semimajor axis was already sufficiently large. 

Figure [3] shows f]{t) in the case where ap = 0.675 and 
Ep = 0.3. The variations due to the planet's orbit are clearly visible 
as rapid fluctuations with amplitude of approximately 0.02. While 
eccentricity is the principal source of fluctuations here, the influ- 
ence of Q Cen A is clearly visible (e.g. at t ~ 350 years), reducing 
the habitable fraction by ~ 0.01. In this example, equatorial lat- 
itudes exhibit temperatures above the boiling point of water - the 
passage of a Cen A increases the thickness of this inhospitable 
band. 

3.3 Binary Transience 

If we now take the previous example and set Bp = 0, then we can 
see the influence of a Cen A begin to dominate (left panel of Fig- 
ure |4l(. The fluctuations due to eccentricity have effectively disap- 
peared, and periastron passage of a Cen A (right panel of Figure|4] 
induces sharp decreases in fj{t) of around 0.025. This clearly shows 
the habitability of the planet being significantly affected by a Cen 
A. This is despite the ratio of mean insolation from the primary and 
secondary, ^ < 0.01. 



4 DISCUSSION 

4.1 Dependence on the Secondary's Orbit 

We have focused specifically on the a Cen system in this paper. It 
is instructive to investigate other putative binary systems, to investi- 
gate the varying strength of binary perturbations. As the parameter 
space of binary star systems is quite large, we restrict ourselves to 

The same value is achieved at ap = 0.69 au, but this solution has no ice 
caps, and an equatorial T > 373 K 



varying the orbital parameters of a Cen A, (a2,e2), and assume 
the entire system is coplanar. Figure [5] shows how a planet's min- 
imum, maximum and global mean temperature varies as a2 is de- 
creased (with 62 held at 0.5179). The planet parameters are fixed at 
ap = 0.9, Cp = 0. 

Taking the true orbital parameters of a Cen A (bottom right 
of figure), we can see that this is indeed a snowball planet, with 
the maximum temperature never exceeding 210 K. A slight per- 
turbation in the mean can be seen at f ~ 350 years due to pas- 
sage of a Cen A through periastron. Decreasing 02 also decreases 
the orbital period; this can be seen in the other panels of Figure 
[S] with the minimum a2 = 5 au (top left). More periastron pas- 
sages induce more temperature fluctuations, but these are unable 
to melt the planet from its snowball state, with the exception of 
02 = 5 au (which produces mean temperatures close to terrestrial 
values). This climate cycle is stable over the entire simulation time. 
In this configuration, the insolation produced by a Cen A consti- 
tutes nearly a third of the total insolation by both stars at periastron, 
and around ten percent of the total insolation at apastron, so it is not 
entirely surprising that the planet can escape the snowball state. 

Whether the outer HZ is extended or not depends quite sensi- 
tively on a2. We can see that 02 = 6 au (top right) fails to melt the 
planet, despite producing temperature fluctuations of order 10 K. 
The insolation due to a Cen A is still approximately a third of the 
total energy budget at periastron, but the reduced frequency of pe- 
riastron passages due to orbital distance means the time-averaged 
insolation is lower. This very small change in the magnitude and 
frequency of the perturbation induced by a Cen A is sufficient to 
completely alter the state of the planet in orbit around a Cen B 
from habitable to snowball. 

In the habitable case shown here, the periastron of a Cen A 
is approximately 2.4 au. Can we produce another habitable planet 
by maintaining this periastron and a Cen A's true semi major 
axis, a2 = 23.23 au? This would correspond (roughly) to an ec- 
centricity of 62 ~ 0.9. As the periastron passage is now quite 
rapid, we should also expect that the climate model will become 
much more sensitive to the phase between the planet's orbit and a 
Cen A's. We investigate the dependence on phase by running the 
same simulation with different initial planetary orbital longitudes 
(0p = 0°, 90°, 180°). 

The results can be seen in Figure[6] We can see that the planet 
remains in a snowball state, but experiences temperature fluctua- 
tions as large as 20 K. The strength of the fluctuations clearly varies 
with the relative position of planet and secondary as periastron is 
reached - each of the three plots shows the same behaviour, but 
with a phase shift caused by the planet's shifted starting position. 
The perturbation in insolation due to a Cen A retains its maximum 
value at insolation of around 30% of the total, but the orbital period 
remains too long for these close passages to be sufficiently frequent 
to melt the planet. 

4.2 Dependence on Ocean Fraction 

We have assumed that the planets have the same surface ocean 
cover as Earth. Planets with lower ocean fractions are subject to 
shorter thermal relaxation ti mescales, and there fore will experience 
greater seasonal variations jSpiegel et al.ll200§) . It might therefore 
be reasonable to assume that drier planets will be more sensitive to 
binary fluctuations than wetter planets. This indeed appears to be 
the case, as is shown in Figure |7] We display habitability fractions 
for the same planetary parameters as Figure|4]with f ocean now 0.1. 
The seasonal variations are substantially higher, with the fluctua- 
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Figure 4. Left:Habitability fraction as a function of time for tlie case wliere = 0.675, ep = 0.0. Riglit: distance between tlie secondary and the planet as a 
function of time for the same planet parameters. 
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Figure 5. Comparing planet temperature as a function of time as the semi major axis of a Cen A is changed (for a planet with ap = 0.9, e = 0). Top left 
shows a2 = 5AU, top right shows a2 = 6AU, bottom left shows a2 = lOAU, and bottom right shows the true a Cen A semi major axis of 23.23 AU. The 
solid line indicates the global mean temperature, the dashed line the minimum temperature, and the dotted line is the maximum temperature. 



tions due to a Cen A also slightly increased (from 0.025 to about 
0.03). 



4.3 Dependence on Obliquity 

The influence of a Cen A adds a second seasonal variation to plan- 
ets around a Cen B, with the polar declination varying in a non- 
trivial fashion over the period of the binary's orbit. How will binary 
oscillations alter under a change in obliquity? We have used the ter- 
restrial value of 23.5° until this point, but there is no reason to dis- 
count other values. In our Solar System, Mars' obliquity appears to 



have varied significantly between 0° and 60° fcaskar & Robuteil 
Il993h . In the case of Earth, the Moon has played an important 
role i n stabilising obliquity fluctuations In eron de Surgv & LaskaJ 
but this is not necessarily a general result for terrestrial plan- 
ets with relatively large moons. 

The influence of a binar y can help lock the planet's obliquity 
into a so-called Cassini state jCorreia et al.ll20lTh . While these are 
generally low obliquity states, high obliquity states can also oc- 
cur (JDobrovolskis 2009). Even without binary influence, numerical 
simulations indicate that the distribution of terrestrial planet obliq- 
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Figure 6. Comparing planet temperature as a function of time witli the eccentricity of a Cen A, boosted to 62 =0.9 (for a planet witli ap = 0.9, e = 0). Top 
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uity is iso tropic, and ther efore primordial obli quities may be large 
jKokubo & Ida 2007; Mi guel & Bruninill2010l) . 

ISpiegel et al.l ( |2009|) show that high obliquity planets experi- 
ence stronger seasonal variations, and can move far from global 
radiative balance even in circular orbits. While they are not neces- 
sarily more prone to snowball transitions, we may expect them to 
be more sensitive to binary fluctuations. 

We investigate two other values of obliquity for the case where 
Op = 0.675 au, = 0, increasing the obliquity to 45° and 90° 
(Figure [8](. The amplitude of the temperature oscillations are un- 
affected, but the temperatures themselves change significantly. In- 
deed, in the case of 45° spin, the planet requires re-classification 
from a binary tran sient to a habitable planet. This change (in line 
with the results o f ISpiegel et al1l2009l for single-star systems) un- 
derlines the difficulty of classifying planets as habitable or other- 
wise, as the parameter space for habitability is non-trivial and high 
in dimension. 



4.4 Dependence on Rotation Rate 

Other than convenience, there are no real reasons to assume that 
planets around a Cen B possess the same diurnal period as Earth. 

We investigate planets rotating with 8, 24 and 72 hour periods 
(with flp = 0.675, Bp = 0). Figure|9]shows the resulting minimum, 
maximum and mean temperatures as a function of time for all three 
cases. The temperature fluctuation induced by a Cen A maintains 
the same amplitude regardless of rotation rate. 

Reducing the rotation period reduces heat transport, causing 
the radiant energy deposited by the primary to be diffused to a 
smaller range of latitudes, increasing the temperature gradient in 
the fast rotating case (top panel of Figure |9}. The slower rotating 
case shows the more efficient latitudinal transport distributing inso- 
lation such that the temperature range from pole to pole is less than 
20K (bottom panel). 
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Figure 8. Comparing planet temperature as a function of time as the obliquity of the planet is changed (for a planet with = 0.675, e = 0). The left hand 
plot shows a planet with with obliquity increased to 45°, and the right hand plot shows a planet with obliquity increased to 90°. The solid line indicates the 
global mean temperature, the dashed line the minimum temperature, and the dotted line is the maximum temperature. 
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4.5 Limitations of the Model 

Finally, we sho uld note the li mitations of the above analyses. As 
with the work o f lSpiegeletai] (l20081, the ID LEBM can only cap- 
ture a restricted range of thermal timescales. While it is sensitive to 
seasonal forcing and (as a result) orbital timescales, the model does 
not exhibit the features of much slower climate processes, which 
possess timescales many times the orbital period of the a Cen sys- 
tem. For example, oceanic circulation is not included in this model, 
which can determine longer-term climate variation. 

More compl ex additions, such as clou ds and a carbonate- 
silicate cycle (e.g. I Williams & KastindfTggTh would provide a fur- 
ther regulating effect on planet temperatures, potentially extend- 
ing the outer edge of the habitable zone (Kasting et al. 1993) on 
long enough timescales. In the particular case of binary systems, 
we might expect that the oscillations induced by the presence of a 
companion would be damped, if the silicate cycle can respond suf- 
ficiently quickly during periastron passage. The timescale on which 
the silicate cycle can be expected to respond is sensitive to the spe- 
cific properties of the planet, i.e. its chemistry, geology and ocean 
circulation systems. For the Earth, it has been estimated that the 
equil ibration timescale for the ca rbonate-silicate cycle is around 0.5 
Myr ^Williams & Kastin3ll997h . A more rapid means of CO2 re- 
lease may be through warming the oceans, as their capacity to hold 
CO2 decreases with increasing temperature. The relevant timescale 
would instead be ocean circulation timescales, which are of the or- 
der 10'' years. Again, these timescales are representative for the 
Earth only, and will not apply in general to terrestrial exoplanets. 

Perhaps the most obvious limitation of these models is their di- 
mension. ID modelling prevents discussion of longitudinal climate 
variations, and forces us to consider diumally averaged insolation. 
This prevents simulation of planets with sl ow rotation rates re lative 
to their orbital motion (such as Venus, see jParish et al.ll20Tlh . The 
contrast between land and ocean is also lost in ID, and latitudi- 
nal variations that occur as a result are not accounted for. Despite 
this, ID LEBMs still capture much of the relevant physics, capa- 
ble of reproducing fiducial Earths wit h temperature profiles very 
similar to real data jSpiegel et al.l2008l) . The inner edge of the hab- 
itable zone is less well-defined than the outer edge - atmospheric 
changes could allow liquid water above 373 K, and the runaway 
greenh ouse effect m ay becoine important at temperatures nearer 
350 K dSpiegel et alj ,2008 and references within). In any case, the 
outer edge is likely to be more interesting from an astrobiological 



standpoint, as current and future instrumentation will be more ca- 
pable of probe spectral features of pl anets at larger semi-major axes 
(see e.g. lKaltenegger & Selsisl201U) . 

We should acknowledge, however, that the addition of a sec- 
ond insolation source into diffusion approximation-based models 
such as this is not entirely understood. However, the perturbations 
induced by secondary insolation are relatively small (typically a 
few percent of the primary insolation, except in extreme cases). We 
therefore argue that the models are appropriate for an initial ex- 
plor atory investigation. 

iThevenin et al.l ( |2002[) estimate the luminosity of a Cen B as 
0.5002 ± 0.016Z/Q, which differs from the main sequence value 
used in this work by around 30%. The habitable zone would be 
somewhat closer to the star as a result, increasing the number of 
planetary orbits per binary orbit. While much of the qualitative 
trends in this paper would remain unchanged, it is clearly important 
that future studies of habitability in this system use observationally 
constrained luminosities. 

We should also note that we neglect the effect of a Cen A's 
gravitational field on the dynamics of a Cen B's planets. While 
the planet's orbital inclination is low enough to avoid Kozai reso- 
nances, and the mean motion re sonances due to a Cen A are absen t 
from the orbits we explore ( Mic htchenko & Porto de Melldl2009h , 
we should still expect short period perturbations to the planet semi- 
major axis of order 0.01 au, more than sufficient to produce sub- 
stantial Milankovitch cycles jSpiegel et al.ll201oh . We have seen 
that the outer edge of the HZ is sensitive to perturbations of this 
size, so future studies must incorporate these perturbations in more 
detail. 



5 CONCLUSIONS 

We have investigated the influence of a Cen A on the habitable 
zone of a Cen B, to test the efficacy of the single-star approxima- 
tion in a binary context. In general, we demonstrate that the single- 
star approximation is roughly correct for calculating the inner and 
outer boundaries of the habitable zone, but fails to capture oscil- 
lations in the planet's climate that occur as a result of a Cen A's 
passage through periastron. 

For all planetary orbits, the presence of a Cen A induces tem- 
perature fluctuations of order a few K. At the habitable zone bound- 
aries, the fraction of habitable surface on such planets can be altered 
by around 3%. The strength of these fluctuations can be increased 
by reducing the planet's ocean fraction, or increasing its obliquity. 

It is reasonable to speculate that if life were to exist on planets 
around a Cen B, that they may develop two circadian rhythms (cf 
iBreus et^l 19951) corresponding to both the length of day around 
the primary, and the period of the secondary's orbit (approx 70 
years). Altering the available habitat by a few percent may also 
influence migration patterns and population evolution. 

While we have demonstrated that the temperature fluctuations 
for planets around a Cen B due to a Cen A are relatively small, 
the consequences of a periodic temperature forcing of a few K to 
long term climate evolution cannot be fully understood from this 
work. To fully appreciate the impact on (for example) ocean cir- 
culation and carbonate-silicate cycles requires further investigation 
with more advanced climate models. 
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